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In correcting satellite doppler data for tropospheric effects, 
we have found that we can infer the precipe table water vapor (PWV) 
at the tracking site. The determination proceeds via a least-squares 
fit of a refraction parameter to each fifteen-minute pass of data. 

For its success, the technique depends on having: 1) an ephem- 

eris for the satellite, 2) an analytic model for the refraction range 
effect that is good to a few centimeters, 3) doppler (range differ- 
ence or cycle count) data with noise level below 10 centimeters, and 
4) a surface pressure/temperature measurement at the tracking site. 

The PWV is a by-product of the computation necessary to correct 
the doppler data for tropospheric effects. We have tried for a num- 
ber of years to reliably isolate the tropospheric refractive effect 
in the doppler shift measurement. It is only recently that we have 
succeeded in doing this. Our rec.r.t success was due to a formulation 
of the refraction integral which minimized the necessity for explicit 
water vapor, temperature and pressure profiles. 

1. THE DOPPLER SHIFT AND THE REFRACTIVE EFFECT 

A common and convenient model of the doppler shift measurement 
contains the range to the satellite at two contiguous times. These 
ranges are corrupted by tropospheric effects which can be written: 


(N.+N ) dp (1) 

d w 

wherein 

N^ is the ’dry" refractivity 
N w is the "wet" refractivity 
The Smith-Weintraub expressions for these are: 
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N d =77.6 P/T 

N = (3.73 x 10 5 °K 2 /mbar ) e/T 2 
wherein 


(e,T,P) are the (water vapor pressure, temperature total atmo- 
spheric pressure) at a point within the troposphere. 
Units consistency requires e and P in millibars and T 
in °K 


Physical rigor (Fermat's principle) insists that the integration path 
be an extremum. Except at low elevation angles, this requirement is 
not important as ray-tracing studies show; we can use the instanta- 
neous straight line connecting the satellite and observer for the 
integration path. If we then change variables from p to h (see 
Fig. 1) , we obtain 



(N ,+N ) 
d w 



( 2 ) 


Since both the wet and dry troposphere extend to heights that are 
small compared to the radius of the earth,- h/Re < h . , = 0.007, we 

replace eq. (2) with ' e 



(N.+N ) dh 
a w 


wherein 


dw 


H 




cos E 

T 


- 1/2 


dw d/Re, 


( 3 ) 


( 4 ) 
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Using explicit models of and N w which, in turn, require models for 
the atmospheric state, we were able to show that eq. (3) and (4) 
could be written 


As = 


10 


-6 


/ (K 


.+11 ) dh 
d w 


l -I c o s, 

1 \1 + 0 . 


— V 

001 ) 


1 - 1/2 


(5) 


with an error which did not exceed a few centimeterr, if we do not use 
elevation angles less than 7.5®. Above 20° the form Is" practically 
exact, i.e., the error is, at most, a few millimeters. 

Equation (5) is a convenient separation (amplitude . geometry) , 
moreover, in the amplitude 


is h 



/* 

J <VV ah 


( 6 ) 


We have absorbed, "buried", the difficult problems of modeling the 
atmospheric state as a function of the surface conditions. We say 
difficult'; the water vapor is not uniformly distributed in the half 
space above the observing site, and modeling has fundamental limita- 
tions. 


We used eq. (5) as a model of the tropospheric range effect and 
associated an amplitude, eq. (6), with each fifteen-minute pass of 
data. We fitted the amplitude without any a priori constraints on the 
value of As^; although we have, with a surface pressure measurement, a 
lower bound" 


(the dry effect = 10 



dh) which is 85-90 percent of the 


combined wet-dry integral. This dry integral is very closely 2.305 
(meters) times the surface pressure in atmospheres. To isolate the 
tropospheric amplitude, we performed a four-parameter fit; the three 
additional parameters absorbed orbit, position and frequency biases. 


2 . RESULTS 


We have used doppler data characterized as follows: 

1. Obtained at a globally-distributed set of tracking sites 
from two of the Transit satellites, 30140 and 30180. 

2. Ore two-day data span was obtained during February 1980 
(days 48-49) and another during July 1982 (days 197-198) . Each span 
contained approximately eighty passes. 



Results of the amplitude fitting for days 197-198 are shown in 
Fig. 2. The lower bound, the straight line, is a consequence of the 
Hopfield theory (see Refs. 2 and 4) . Data points falling below this 
straight line are "wild" (invalid). There are four of these. The 
distance from each of the valid data points to the straight line is 
the wet term 


10^ f N dh 

J w 

We then computed the mean value (for the valid data) at each of the 
sites; mean over two days for each of the two data sets. Four sta- 
tions were represented in both data sets: Johannesburg, South Africa; 

Las Cruces, New Mexico; San Jose Dos Campos, Brazil; and Herndon, 
Virginia. The data at these four sites agree with the (usually true) 
finding that the atmosphere contains more water vapor in the summer 
than it does in the winter. Fig. 3. A tabulation of the results shown 
in Fig. 3 is given in Table 1. From Fig. 3 we see that synoptically 
the wet term is about 10 percent of the dry. 

3. PRECIPITABLE WATER VAPOR 

The 'wet' range integral for a vertical path can be written, 
using the Smith-Weintraub expression for the refractivity : 


As 

w 


0.373 



% dh 
T 


(7) 


and assuming the water vapor obeys the perfect gas lav;: 


e = p RT 
e 


(8) 


we have 


As 

w 


0.373 




(9) 


and, in turn, a linear decrease of temperature with height facilitates 
an expansion in terms of the moments of the water vapor distribution. 


As w = 


0.373 


Rp 


w 


kf 


p dh + 


v 


■J 


h p^ dh + 


( 10 ) 


alternately, 


331 







ORIR(n AL p age , 3 

OF POOR QUALITY 



nzvua— ■ 




Table 1 

Troposphere Fitting Results 
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D w is the density of (liquid) water. Tq is the surface temperature 
( °K) and L is the (absolute value of the) temperature lapse rate 
(- 6.5°/km), h wv is the scale height of the water vapor in the atmo- 
sphere (- 2.2 km). Each of the teriris in the expansion, eq. (10), has 
units of height. T H e first term, the dominant one, is the "precipi- 
table water vapor". The dimensionless term in front is about 6.0. 

It is precisely 6.0 for Tq = 287°K. The second term (the first moment 
in the expansion) is approximately 5 percent of the precipitable water 
vapor (Ref. 3) . The second moment is a factor of ten less than the 
first. Consequently, we can, with surface temperature and pressure 
measurements, interpret our results for the (wet) range effect in 
terms of precipitable water vapor. 

4 . CONCLUDING REMARKS 

In keeping with the spirit of this meeting, the significant find- 
ings from this study are: 

1. It is possible to determine the total atmospheric water vapor 
using microwave doppler measurements. This measurement differs from 
the water vapor measurements made with radiosondes or with microwave 
radiometers. The latter two produce water vapor estimates along a 
profile through the atmosphere. The doppler technique produces a 
different sort of estimate; it is a 'bulk' estimate over that fraction 
of the troposphere swept out by the satellite-observer vector. 
Typically a region extending out to 150 km, on the satellite side of 
the observing site, will be sampled (see Ref. 1) . Each of these kinds 
of measurements have their uses. 

It would be interesting and useful to study the correlation 
between the different kinds of measurements. The data currently 
exists to do this. 

2. It is possible to derive a model of the tropospheric range 
effect which is accurate to a few centimeters. The 'thinness' of the 
troposphere (compared with the radius of the earth) is an essential 
fact in developing the model. By model we mean one in which there is 
a clear separation of geometry on the one hand and an atmospheric- 
state-dependent amplitude on the other. 

3 . The tropospheric range amplitude is (demonstrably) separable 
from the other errors affecting the doppler shift. The accuracy with 
which the tropospheric range effect is isolated depends on the 
accuracy of the modeling and the noise level of the data. The exist- 
ing data noise is about 10 cm (in range difference) , and we are cur- 
rently unable to fully exploit the accuracy implicit in the theory. 
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It is easily v’ithin the current state-of-the-art to build 
satellites and ground equipment which will reduce the noise level of 
the data. Equally important for these purposes would be raising the 
frequency pair from 150/400 MHz to (say) 400/1200 MHz to reduce the 
higher order ionospheric effects in the data. 
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